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' The structure and dynamics of atomic oxygen adsorbed on Ag(410) and Ag(210) surfaces has been 

, investigated using density functional theory. Our results show that the adsorption conhguration in 

\^ • which O adatoms decorate the upper side of the (110) steps forming O-Ag-O rows is particularly 

' stable for both surfaces. On Ag(210), this arrangement is more stable than other configurations 

at all the investigated coverages. On Ag(410), adsorption on the terrace and at the step edge are 
almost degenerate, the former being slightly preferred at low coverage while the latter is stabilized by 
' increasing the coverage. These findings are substantiated by a comparison between the vibrational 

modes, calculated within density-functional perturbation theory, and the HREEL spectrum which 
has been recently measured in these systems. 

O ' PACS numbers; 68.43.-h 

^ I. INTRODUCTION 

The interaetion of oxygen with silver surfaces has been intensively studied because of the key role of silver in 
' many important industrial oxidation reactions such as, for example, partial oxidation of methanol to formaldehyde 
or ethylene epoxidationi These reactions occur very efhciently when Ag powders are employed as a catalyst. Since 
the particles which compose these powders are affected by many surface defects of different nature and with different 
concentrations, their chemistry can be extremely different from that of perfect low-Miller-index surfaces. Indeed 
Q • the presence of various inequivalent adsorption sites suggests that atomic and molecular oxygen can adsorb forming 
^ ' species of different chemical nature. The microscopic characterization of such adsorbed species is a fundamental issue 
to understand the catalytic activity of Ag powders. 
I ' A large amount of experimental and theoretical results clearly indicates that steps, kinks, and other surface imper- 
^ , fections play a crucial role in surface chemistry. For example, the study of many systems such as N and O adsorbed 
I ' on Ru(0001)(^ O on Pt(lll)r' or C2II4 on Ag(lOO)^ has shown that monoatomic steps can bind adsorbates more 
t — " strongly than terrace sites. Moreover it has been shown that this kind of defects can completely determine the kinet- 
I ics of dissociation in systems such as, for example, N2 and NO on Ru(0001)f^ O2 on Pt(lll)fi N2 on Fe(llO) and 
^ ; Fe/Ru(0001).^ 

The role of steps in the dynamics of 02/Ag interaction has been recently investigated by Rocca and coworkers2ii2iii 
who have studied the Ag(410) and Ag(210) surfaces. These surfaces are vicinal of Ag(lOO) and are characterized by 
open (110) steps and three-atom- and one-atom-row wide (100) terraces, respectively. Using a supersonic molecular 
' beam to dose O2 onto Ag surfaces at selected angles of incidence and different kinetic energies, and characterizing the 
adsorbed species by vibrational spectroscopy (HREEL), these authors have investigated the various oxygen species 
" ■ which form on these surfaces. It was concluded that O2 dissociation occurs mainly at the steps and that atomic 
■ O ' oxygen adsorbs at different surface sites. In particular, HREEL spectra show two peaks in the frequency region of the 
C . 0/Ag stretch modes, at about 32 and 40 meV. These frequencies are quite similar to those observed in O/Ag(100) 
^ ' (33 meV)i^ and on the added-row reconstructed O/Ag(110) surface (40 meV)>i^ It was thus proposed to assign the 
low energy peak to adatoms occupying (100) terrace sites and the one at higher energy to atomic oxygen adsorbed at 
the step edges forming O-Ag-O chains similar to the added rows found in reconstructed O/Ag(110) surfacesi ^^d^ In 
O/Ag(210) an additional peak around 56 meV was found, and it was proposed that it is due to the vibration of an 
5^ • oxygen atom occupying a subsurface site. 

' Motivated by these results we have performed ah initio calculations aimed at understanding the role of the steps in 
the adsorption of atomic oxygen on the Ag(410) and Ag(210) surfaces. In the present work we limit our investigation 
to on-surface atomic oxygen, subsurface and molecular adsorption will be addressed in separate papers. Our results 
show that on these vicinal surfaces the adsorption configm-ation in which the adatoms decorate the steps in a (1x1) 
geometry is particularly stable. The oxygen adsorption on the two surfaces presents however significant differences 
related to the availability of adsorption sites on terraces. Indeed, on Ag(410), terrace and step edge sites have 
comparable chemisorption energies and, in particular, the hollow sites on terraces are slightly more favored than 
the step edge sites when the adatoms are far apart from each other. At higher coverage, instead, adatoms slightly 
prefer to decorate the steps. On Ag(210); the absence of favorable terrace sites makes the step decoration much 
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more stable than any other adsorption configurations. The formation of these 0-Ag-O chains at step edges will be 
shown to significantly affect both the geometrical structure and the electronic properties of the surface. As we have 
already mentioned, the occurrence of adsorption on both the terraces and steps of vicinal Ag(lOO) surfaces was also 
proposed based on an analysis of the HREEL spectra, recently measured for these systems Aifliii Our calculation 
of the vibrational properties of these systems — performed within density-functional perturbation theory — provides a 
further support to this analysis. 

The paper is organized as follows. In Section ^] we describe the theoretical approach and computational details of 
our work. Section IIIII contains our results for the energetics, the geometrical and electronic structures, and for the 
vibrational properties of the systems we have investigated. The last Section is devoted to our conclusions. 

II. COMPUTATIONAL FRAMEWORK 

All the calculations have been performed within density functional theory (DFT) using the generalized gradi- 
ent approximation (GGA) with the Perdew-Burkc-Ernzerhof (PBE) exchange-correlation functionalji^ Vibrational 
properties have been calculated using density-functional perturbation theory (DFPT)iiii^ii2i2£ We have used the 
pseudopotential method with ultra-soft pseudopotentialsSi and plane-wave basis sets up to a kinetic-energy cutoff of 
27 Ry (216 Ry for the charge-density cutoff). Details about the Ag and O pseudopotentials are reported in Ref. 12^. 
Brillouin zone (BZ) integrations have been performed with the Gaussian-spreading special-point techniquo^iS^ with 
a smearing parameter of 0.03 Ry. All the calculations have been performed using the PWscf and PHONON packagesjS^ 
while molecular graphics has been generated with the XCRYSDEN package]^ 

Surfaces are modeled with periodic super-cells. For the Ag(410) and Ag(210) surfaces we have used slabs of 20 
(410) layers and 14 (210) layers, respectively. For the Ag(lOO) surface and the (2x1) added-row reconstructed Ag(llO) 
surface [Ag(110)p(2x 1)0] we have used slabs of 7 layers. Adjacent slabs are separated by a vacuum region of at least 
16 a.u. and O atoms are adsorbed on both sides of the slabs. All the structures have been fully relaxed until 
the Hellmann-Feynman forces were smaller than 10~^ Ry/a.u. Atomic oxygen adsorption on Ag(410) and Ag(210) 
has been modeled by (1x1) and (2x1) surface super-cells. For the Ag(410) surface the BZ integrations have been 
performed using a (8x4x1) and a (4x4x1) uniform shifted fc-mesh2i for the (1x1) and (2x1) structures, respectively, 
while for Ag(210) (8x7x1) and (4x7x1) meshes have been used for (1x1) and (2x1) structures, respectively. For 
the c(2x2) cell on Ag(lOO) we have used a (9 x 9 x 1) mesh, while for the Ag(110)p(2xl)O surface we have used a 
(6 X 9 X 1) mesh. Spin polarization effects have been neglected after checking that the interaction between O and Ag 
substrate results in a negligible spin moment on the O atom. 

The chemisorption energies, E'chem, are referred to the clean Ag(nlO), n = 4 or 2, surface and the isolated oxygen 
molecule: 

^chcm = {Eo/Ks - - NoiEoj2))/No, (1) 

where the total energy of the adsorbate-substrate system, of the clean surface, of the isolated O2 molecule, and the 
number of adsorbed oxygen atoms are represented by i?o/Agi EAg, Eq^ and Nq, respectively. With this definition, 
stable adsorbates have negative chemisorption energies. Considering the errors in the structural relaxation and the 
energy convergence with respect to the plane wave cutoff energy and the k-point density, we estimate an overall 
numerical accuracy for the chemisorption energies of about 30 meV. 

The PBE GGA functional used in the present work slightly overestimates the binding energy of the O2 molecule: 
5.50 eV, as calculated using spin-polarized GGA with a 15 A-wide cubic box, while the experimental value is 5.23 
eV^ Molecular overbinding is a common drawback of current implementations of DFT. This fact, of course, does 
not affect the conclusions of the present work which mainly concerns the relative stability of different adsorption 
configurations. 

Vibrational frequencies are calculated by diagonalizing the q = dynamical matrices of large slabs containing more 
that 50 atomic layers. The common procedure for calculating the dynamical matrices of such large slabs is to patch 
the dynamical matrices obtained for much smaller slabs (in fact, the same which were used to determine the structural 
properties, as described above) with the interatomic force constants calculated for the bulk metaliSS The dynamical 
matrices of the Ag bulk have been calculated on a 4x4x4 grid of q-points in the BZ and the interatomic force 
constants have been obtained by a Fourier transform^ For the thin slabs, we have computed only the parts of the 
dynamical matrices which refer to the displacements of the adsorbates and of the uppermost Ag layers. Diagonalizing 
the dynamical matrix of the extended slab we have calculated the vibrational modes of the systems and we have 
identified the adsorbate-substrate modes from the displacements eigenvectors. The estimated numerical accuracy of 
the vibrational frequencies so obtained is of about 1 meV. 
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III. RESULTS 
A. Energetics 

In Fig. n]we display the four adsorption sites that we have considered in the present work: site A with O just above 
the step and coordinated with three Ag surface atoms; sites Ti and T2 with the adatom in the hohow sites on the 
(100) terrace; site B where O adatom lies in the hollow site just below the step. Note that no terrace Ti and T2 sites 
are available on the Ag(210). Top and bridge sites are not considered here since we have found that they are much 
less stable than hollow sites, as on Ag(lOO) surfacei^ 

1. 0/AgUlO) 

Let us first consider the adsorption of atomic oxygen on the Ag(410) surface. In TablcQlwe present the chemisorption 
energies corresponding to different adsorption configurations (see Fig. |2l. At low coverage (0 = 1/8 ML)^ where 
the adatoms stay far apart from each other, the hollow sites on terraces (Ti and T2) and the step-edge site (A) are 
almost degenerate in energy, with terrace sites being only slightly more favored than site A, while site B is much 
less stable. In these configurations the distances between the O atoms are quite large (8.32 A) and this suggests a 
negligible interaction between the adatoms, as confirmed by test calculations done at lower coverages in which the 
0—0 distance is increased. At higher coverage (O = 1/4 ML) we have examined various possible geometries which 
we indicate with the notation 'S1-S2', where the S's stand for two near A, B, Ti and T2 sites being occupied by two 
O adatoms in a (2x1) super-cell. The Si-Si arrangements correspond to a (1x1) periodicity where O adatoms form 
rows parallel to the steps. Those configurations in which the separation between the adatoms is larger than 5 A are 
not considered since we expect that their adsorption energies should be very similar to those at low coverage. Indeed, 
as we will show, already for 0-0 distances as small as 4 A the interaction between the adsorbates has small effects on 
the chemisorption energies. Our results show that the Si-Si configurations with both the adatoms on the terrace or 
above the step and the A-T2 arrangement are the most stable adsorption geometries. As in the case of low coverage 
the differences in chemisorption energy between these sites are quite small, but in this case we find that oxygen 
adatoms slightly prefer to decorate the step edges. One expects that occupying adjacent sites should be unfavorable 
because of the electrostatic repulsion between negatively charged oxygen adatoms. Indeed, for example, Feibelman 
and co-workers^ have shown that O adatoms on Pt(211) vicinal surface decorate the steps in a (2x1) adsorption 
geometry while the (1x1) structure — where the distance between adjacent oxygens is about 2.8 A — is significantly 
less stable. Our results indicate that this effect on Ag(410) surface is important only when nearest neighbor sites 
are occupied, i.e. when the distance between the adatoms is less than 4 A (as in A-Ti, T1-T2, T2-B and A-B, 
configurations), while it is weak and uneffcctive in the other geometries. Compare for example the A-T2 and A-Ti 
configurations which consist of one adatom sitting above the step while the other sits at a terrace site. In the A-T2 
arrangement the distance between the O atoms is 4.23 A and their chemisorption energy (—0.80 eV) is quite close 
to the values found for oxygen adsorption at sites A and T2 at low coverage (—0.75 eV and —0.83 eV), indicating a 
small adsorbate-adsorbate interaction. In the A~Ti configuration, instead, the 0-0 distance is only 3.26 A and the 
adsorption energy of the adatoms (—0.59 cV) is strongly reduced with respect to when oxygen is adsorbed at sites A 
and Ti at low coverage (—0.75 eV and —0.80 eV). Moreover, it is interesting to observe that in the A-Ti geometry 
the electrostatic repulsion strongly affects also the locations of the adatoms. Indeed we find that the adsorbates are 
slightly displaced from the hollow positions and that their distance, 3.26 A, is considerably larger than the separation 
between two nearest hollow sites of the clean surface (about 2.9 A). This is a clear indication that the two adatoms 
repel each other. 

2. O/Ag(210) 

Similar calculations have been performed also for the Ag(210) surface. The two surfaces differ in the width of 
the terrace (see Fig. No terrace hollow sites are available on the Ag(210) surface. In Table UTI we report the 
chemisorption energies for O adsorbed in various configurations on this surface. These results display the same trend 
as obtained for Ag(410), although the chemisorption energies are somewhat smaller here. Due to the lack of stable 
adsorption sites on the terraces, the formation of 0-Ag-O rows at the step edges (i.e. the A-A configuration) is 
by far energetically favored on this surface. The second most stable configuration at higher coverage (8 = 1/2 ML) 
is therefore A-B, which is 0.24 eV less bound than A-A. In contrast, on Ag(410) several other configurations that 
involve terrace sites are more stable than A-B and much closer in energy to the A-A configuration. Comparing the 
chemisorption energies of various configurations, irrespective of the O coverage, we find that the difference between 
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the most favored A-A configuration and the second most stable one is 0.03 and 0.12 eV for Ag(410) and Ag(210) 
surfaces, respectively. 

3. The stability of 0-Ag-O rows 

One of the most peculiar results of our calculations is the increase of the stability of oxygen adatoms when they 
decorate the step edge in a (1x1) geometry. We believe that the mechanism responsible for this behavior could 
be related to the arrangement of O adatoms to form 0-Ag-O chains at the upper side of the (110) steps. Indeed, 
the A-A adsorption configuration is similar to that occurring in the reconstruction of the Ag(llO) surface upon 
oxygen adsorption, in which O atoms align in added-rows along the [001] direction, thus forming stable 0-Ag-O 
structuresi^^ (the Ag(110)p(2x 1)0 surface is shown in Fig. ISJ. As we have pointed out in the introduction, this 
similarity between the two structures was first suggested from a comparison of the vibrational properties of the two 
surfaces ii£ 

Before analyzing in detail the structural and electronic features related to the formation of these chains, it is 
interesting to consider the effects of the substrate relaxation on the chemisorption energies. For this reason we have 
repeated some calculations for O/Ag(410) keeping the substrate fixed in the optimized clean surface geometry and 
allowing the adatoms to relax (see values in parentheses in TableP). Our results show that at low coverage the T2 site 
is more stable than A site even if substrate relaxation is not considered. On the other hand, in the (1x1) periodicity 
O adatoms in A-A and T2-T2 configurations have the same chemisorption energy on a fixed substrate. This result 
indicates that the enhanced stability of atomic oxygen on the step in the (1x1) geometry is mainly determined by 
the energy gain due to surface relaxation allowed by the presence of the steps. 

B. Geometrical structure 

In order to illustrate the characteristic features of the 0-Ag-O chain formation at the step edges, we discuss now 
the effects of oxygen adsorption on the geometry of the surface. To facilitate the discussion, we first define a few 
labels for silver atoms — Ag^ , Ag^, and Ag^ — as in Fig. 0pi 

1. 0/AgUlO) 

In table Hill we report the optimized structural parameters describing the oxygen adsorption on Ag(410) at G = 
1/8 ML and at 8 = 1/4 ML for the Si-Si configurations. First of all it is interesting to observe the dependence of 
the bond lengths between the adsorbate and the nearest metal atoms on the number of neighboring substrate atoms. 
As predicted by effective medium theory^^^ the larger the number of nearest metal atoms the longer the bond length. 
Indeed the results in Table UTTl (observe in particular do-Ag" ^-i^d do-Ag^) show that the three- fold coordinated 
step-edge is the site in which O is closest to the surrounding silver atoms; O in the two four-fold coordinated terrace 
sites has similar structural features, while O — when sitting in the site just below the step, close to four atoms of the 
terrace and to two atoms of the step — is the furthest from neighboring metal atoms. 

We now focus on the structural features of oxygen adsorption at sites A and T2. In Fig. [Slwe compare the structure 
of the first three layers of the substrate for the A and A-A geometries. At low coverage the main substrate distortions 
due to the presence of the adatom are a shift of 0.18 A of the nearest Ag terrace atom towards the center of the 
terrace (see c^Agi^-AgS' Table llV|l and an expansion of 0.16 A of the distance between the two Ag atoms on the step 
near oxygen (see dAgR-AgR in Table ITV)l . The result is that the adatom sits in a site equidistant (do-Ag = 2.21 A) 
from the three Ag surface atom surrounding it. In the A-A configuration the distortion has different features. In this 
case symmetry does not allow the two Ag step atoms near oxygen to be displaced parallel to the step, so that these 
atoms move only outward. This is indicated by a marked expansion of the bond lengths between the Ag step atom 
and the Ag atoms on terrace (by 0.25 A) and below the adatom (by 0.16 A). In this case the oxygen is closer to the 
Ag step atoms (do-Ag" = 2.10 A) than to the nearest Ag atom on terrace (c^o-AgS' — 2.24 A). Note also that the 
angle, a, between the bonds of oxygen with the two Ag atoms of the step — which directly measures the alignment 
of the Ag-O-Ag structures — is larger at high than at low coverage (164° versus 156°). Thus, in the A-A geometry 
oxygen and silver atoms on the step tend to align forming an almost straight 0-Ag-O chain which is more stable than 
other adsorption configurations. As we have already pointed out, the 0-Ag-O chains at the step edges have a similar 
structure as the added rows in the reconstructed O/Ag(110) surface, so it is interesting to compare the geometries of 
the two systems. In Table HTTl we report some predicted structural quantities for the Ag(110)p(2x 1)0 surface. Our 
results are in good agreement with experimental, and previous DFT calculations^ As we can observe, the alignment 
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of O and Ag atoms in the added rows on Ag(llO) is more pronounced than in the A-A configuration (compare, for 
example, the angle a in the two geometries), but the similarity between the two structures is clear. 

Fig. Elshows the structures of T2 and T2-T2 configurations. In these systems the height of the adatom with respect 
to the terrace is larger than when the oxygen is chemisorbed on the step. Also in this case in the (1x1) geometry the 
silver atoms in the row of oxygens move up towards adatoms, but here the alignment between O and Ag atoms is less 
pronounced than when oxygen atoms are on the step. 

2. O/Ag(210) 

In Table Ivlwc report the optimized structural parameters describing the oxygen adsorption on Ag(210) at O ~ 1/4 
ML and at = 1/2 ML for the Si-Si configurations, while Fig. Elshows the side views of the corresponding optimized 
O/Ag(210) structures. These data and plots reveal the same trend as for Ag(410) surface. The most apparent feature 
is that the height of the O adatoms above the surface is much smaller for A and A~A sites than for the B and B-B 
sites. Note also the difference between the pattern of the substrate reconstruction for A (low coverage) and A-A (high 
coverage) configurations. Like on the Ag(410) surface, the adsorption of O adatoms induces an outward relaxation of 
step Ag atoms for both the A and A-A configurations. In the case of the A configuration this relaxation is not large 
enough to compensate the inward relaxation of the clean substrate (see the inset of top- left panel of Fig. [7|), therefore 
Ag atoms at the step appear to be pushed somewhat inward, like in the case of the clean Ag(210) substrate. On the 
contrary, this relaxation is quite significant in the A-A configuration, resulting in a marked expansion of the bond 
length between the Ag step and terrace atoms (Ac?^gR_^gS' = 0.20 A), and between the Ag step atom and Ag atom 

located below the O adatom (AdAgR-AgB = 0.19 A). Note that the step Ag atoms appear to be pushed outwards. 
Contrary to the A-A geometry, in the A configuration the step Ag atoms also relax laterally along the step-edge 
direction away from O adatom — the distance between the two step Ag atoms that are bound to the same O adatom 
being increased by 0.17 A. 

C. Electronic structure 

In order to better characterize the different behavior of oxygen on the step in A and A-A configurations we have 
analyzed the corresponding electronic structure. Here we present the analysis of the electronic structure of O/Ag(410) 
only — the conclusions for O/Ag(210) being similar. 

The interaction between oxygen and silver removes electronic charge from the silver atoms neighboring the O 
adatom. The electron deficit regions around the silver atoms are shown in Fig. |Hl(blue isosurfaces) . It is interesting 
to observe that while at low coverage the deficit regions around the three Ag atoms are very similar, thus indicating 
a similar donation of charge, in the A-A configuration the flow of charge to oxygen adatoms is mainly due to silver 
atoms on the step. This suggests that the 0-Ag-O chains that form in the (1x1) geometry have the structure of 
electrostatically stable • ■ • H 1 — • ■ • strings. 

The analysis of the density of states reveals that the bond between oxygen and silver is not purely ionic. In Fig. |5| 
we display the density of states projected (PDOS) onto the oxygen adatom and onto the nearest silver atoms for both 
the A and A-A geometries. The PDOS on O atoms shows that the interaction of the O 2p state with the Ad metal 
band results in the formation of two regions of high density of states, which correspond to bonding and antibonding 
states. The bonding states are located at the bottom of the 4c? Ag band while the antibonding ones are mainly below 
the Fermi level. Note that the O 2s level lies below the bottom of the valence band (at about —17 eV with respect 
to the Fermi level) and it is not shown in the figures. The main difference between the PDOS in the two geometries 
is that at high coverage additional features appear 6-7 eV below the Fermi level. The figure clearly shows that these 
states are due the hybridization of oxygen 2p orbitals and Ad orbitals of the silver atoms on the step (Ag^). By 
inspection of the components of these orbitals we find that the mixing is mainly between O 2px orbital parallel to 
the step-edge and Ag Ad^2_y2 orbital parallel to the terrace plane. This hybridization is very efficient because of the 
good alignment between O and Ag atoms on the step. The spatial distribution of these states is illustrated in Fig. 1101 
where we display the integrated local density of states (ILDOS)^ in the energy window (—7.0, —6.0) eV comprising 
the two small peaks in the PDOS. The figure shows quite clearly the bond between oxygen and silver atoms on the 
step. It is worth noting that in the other Si-Si configurations there are no states in this energy region even if we find 
an increase of the density of states around —6.0 eV, especially in the T2-T2 geometry. 

The study of the electronic structure of the added- row reconstructed Ag(llO) surface confirms that the presence of 
bonding states 6-7 eV below the Fermi level is related to the formation of aligned O-Ag-O-Ag structures. Indeed 
we find similar spectral features also in the DOS of Ag(110)p(2x 1)0 and their spatial distribution closely resembles 
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that shown in Fig. ^| Note that the results of experimental studiesi^*^ show a quite weak oxygen-induced feature 
in this energy range. 

D. Oxygen vibrational modes 

In Table lVlll we present the frequencies of the adsorbate modes and the directions of the corresponding displacement 
eigenvectors (the angles are defined in Fig. Ill() for different 0/Ag systems. Before analyzing the vibrational modes 
of oxygen on Ag(410) and Ag(210); it is useful to discuss the vibrational properties of the O/Ag(100) system with 
oxygen adsorbed in hollow sites and of the added-row reconstructed Ag(llO) surface. Data for the Ag(100)c(2x2)O 
structure — which had already been calculated in Ref. |^in the local density approximation (LDA) — have been recal- 
culated within the present GGA approach, giving quite similar results (not unusually, GGA frequencies result to be 
slightly red-shifted with respect to LDA). The mi mode perpendicular to the surface vibrates at 30 meV in agree- 
ment with the peak at 30-32 meV observed in the HREEL spectra of the non-reconstructed O/Ag(100) surfacei^ 
The modes TO2 and 7713 at 50 meV are two in-plane degenerate modes. For the Ag(110)p(2x 1)0 structure we find a 
vibrational mode, 7712, perpendicular to the surface at 38 meV. This value has to be compared with the feature around 
41 mcV observed in HREEL spectra>i2i^ Besides this dipole-active mode two other adsorbate modes exist, mi at 28 
meV and at 73 meV, in which the adatom vibrates parallel to the surface plane. The polarization of the softer 
mode is perpendicular to the added row while the other one is parallel to it. It is interesting to observe that the ma 
mode is even harder than the in-plane modes (m2 and m^) of Ag(100)c(2x2)O because the oxygen atom sits closer to 
the nearest silver atoms in the added row reconstructed Ag(llO) surface (Az = 0.10 A) than in the Ag(100)c(2x2)O 
structure (Az = 0.71 A)2^. 

We now focus on the vibrational properties of the A-A adsorption geometry on Ag(410) and Ag(210). The mode 
at higher frequency (ma at 68 meV and 71 meV for Ag(410) and Ag(210) respectively) in which the adatom vibrates 
parallel to the Ag-O-Ag chain at the step is very similar to the 1713 mode of the Ag(110)p(2x 1)0 structure. The 
softer modes are due to vibrations of the adsorbate in a plane perpendicular to the step edge (these two modes for 
the Ag(410) surface are shown in the upper panel of Fig. I12|) . The mi mode at 24 meV is nearly parallel to the (110) 
plane of the step while the m2 mode at 37 meV is almost perpendicular to it. These modes resemble very closely the 
mi and m2 modes of the Ag(110)p(2x 1)0 surface. The displacement pattern of these two modes suggests that they 
could both be dipole active. In particular the mode m2 could be responsible of the spectral feature around 40 meV 
observed in the HREEL spectra of these surfaces. 

For the Ag(210) surface we have also investigated the A configuration to understand the effect of the Ag-0~Ag chain 
formation at the step on the vibrational properties. In this case both the frequency and the displacement eigenvectors 
of the softer modes, mi and m2, change very little with respect to the A-A arrangement. The ms mode is polarized 
parallel to the step and is significantly softer than in the A-A structure (46 versus 71 meV) mainly because in this 
case the height of the adatom with respect to the step Ag atoms is larger than in the A-A geometry. 

Since the hollow terrace sites on Ag(410) compete in energy with the step edge sites, we have also investigated 
the vibrational properties of a very stable adsorption configuration which involves terrace sites, namely the Ti-Ti 
geometry. In the lower panel of Fig. E|we show the mi and m2 modes in the plane perpendicular to the step. The 
mi mode at 29 meV is due to the vibration of O nearly perpendicular to the (100) terrace plane while the m2 mode 
at 32 meV is due to O vibrating almost parallel to this plane. The harder mode, m^, is at 58 meV and it corresponds 
to the oxygen vibration parallel to the step edge. Our calculations show clearly that the mi mode is very similar to 
the dipole active mi mode of the c(2x2) O adsorption geometry on Ag(lOO). It seems thus reasonable to expect that 
the spectral feature around 32 meV observed in HREEL spectra on Ag(410) could be due to this stretching mode. 

Finally we have calculated the vibrational modes for the B-B adsorption geometry on Ag(210). An adatom at site 
B is significantly less stable than at the step edge, but it is interesting to investigate its vibrational properties since 
this is the only terrace site available on Ag(210) and the HREEL spectra of this surface show, in addition to the 40 
meV peak, other spectral features around 30 meV and 56 meV. As it is shown in Table IVIII the higher frequency 
mode, m3 at 37 meV, is parallel to the step edge, while the other two modes are in the plane perpendicular to it. In 
particular the mi mode at 27 meV is nearly perpendicular the (110) step plane, while the other, m2 at 34 meV, is 
almost parallel to it. 

Two important results emerge from this analysis. The mi mode in the Ti-Ti configuration of O/Ag(410) at 29 meV 
lies very close to the mi mode of O/Ag(100) (30 meV), in correspondence with the spectral feature around 32 mcV 
observed in HREEL spectra on Ag(410)j'^i^°i" The m2 modes of the A-A configuration in O/Ag(210) and O/Ag(410) 
(both at 37 mcV) lie close to the m2 frequency in the added-row O/Ag(110) surface [Ag(110)p(2x 1)0] at 38 meV, in 
correspondence with the spectral feature around 40 meV observed in the HREEL spectra of O/Ag(410) and O/Ag(210) 
surfaces iSiiSiii We believe that this correspondence between infrared-active vibrational modes, calculated for specific 
configurations, and the features found in the experimental HREEL spectra supports the proposed assignment of those 
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spectral features to atomic species adsorbed on terraces and at step edges. Problems still remain in the interpretation 
of the HREEL spectra of O/Ag(210). On one hand, the feature experimentally found at 32 meV could be due to the 
vibration of oxygen occupying the only available terrace sites (site B), but this assignment remains uncertain since 
the adsorption on these sites is energetically unfavored. On the other hand, no vibration is predicted near 56 meV, 
where a feature in the experimental spectrum is also found. This feature could be the signature of other O species, 
possibly adsorbed subsurface, which we have not investigated in the present work. For sure, these issues call for 
further theoretical and, possibly, experimental work. 



IV. CONCLUSIONS 



In this paper we have investigated the adsorption of oxygen on Ag(410) and Ag(210) surfaces. Our results show 
that the step edge decoration is particularly stable on both surfaces. While on Ag(410) adsorption on terrace and at 
the step edge compete in energy and no clear preference is predictable, on Ag(210) the step decoration is significantly 
favored with respect to other adsorption configurations. The formation of this stable O-Ag-0 rows at the step edges 
strongly affects both the structural and the electronic properties of the surface, and leaves a clear fingerprint in the 
HREEL spectrum of these systems at 40 mcV. 

The situation is less clear for the HREEL peak at 32 mcV. On Ag(410) it is reasonable to assign this spectral feature 
to adatoms occupying terrace sites. On Ag(210) instead this 0-Ag stretch mode could be due to O atoms adsorbed 
on terraces at the step foot, which however are predicted to be metastable adsorption sites. Among the adsorption 
configurations that we have investigated, no candidate was found to support the vibrational mode experimentally 
found at 56 meV. This suggests that other O species, possibly adsorbed subsurface, would be responsible for this 
spectral feature. 
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TABLE I: Chemisorption energies, -Bchcm, and 0-0 distances, do-o, for various chemisorption sites on Ag(410) surface. Values 
in parentheses are obtained keeping the substrate fixed. 



e (ML) 


Configuration 


do-o (A) 


-Bchcm (eV) 




A 


8.32 


-0.75 (-0.66) 


1/8 


Ti 


8.32 


-0.80 




T2 


8.32 


-0.83 (-0.76) 




B 


8.32 


-0.60 




A-A 


4.16 


-0.86 (-0.70) 




Ti-Ti 


4.16 


-0.78 




T2-T2 


4.16 


-0.79 (-0.69) 




B-B 


4.16 


-0.54 


1/4 


A-T2 


4.23 


-0.80 




Ti-B 


4.22 


-0.63 




A-Ti 


3.26 


-0.59 




T1-T2 


3.22 


-0.66 




T2-B 


3.23 


-0.54 




A-B 


3.30 


-0.56 



TABLE II: Chemisorption energies and 0-0 distances for various chemisorption sites on Ag(210) surface. The labels have the 
same meaning as in Tab. Q] 



e (ML) 


Configuration 


do-o (A) 


Schom (cV) 


1/4 


A 


5.10 


-0.68 




B 


5.10 


-0.49 




A-A 


4.16 


-0.80 


1/2 


B-B 


4.16 


-0.42 




A-B 


3.62 


-0.56 



TABLE III: Some key quantities describing the O chemisorption sites on Ag(410) in various geometries. do_AgR is bond length 
with the nearest silver atom, Ag^ see Fig. |1| lying in the row of oxygen adatoms. do_AgS and d^^^^si indicate the bond 
length with the other two oxygen nearest surface silver atoms, Ag^ and Ag^ , while and do_AgB is the bond length with the 
silver atom just below oxygen, Ag^. Az is the height of oxygen with respect to Ag^ atoms and a is the bond angle between 
the bonds with the nearest Ag^ atoms (a — 180° indicates perfect alignment). 



e (ML) 


Configuration 


do-AgR (A) 


do-AgS (A) 


"'o-AgS' (A) 


do_AgB (A) 


Az (A) 


a n 




A 


2.21 




2.21 


2.43 


0.46 


156 


1/8 


Ti 


2.27 


2.26 


2.29 


2.97 


0.80 


138 




T2 


2.28 


2.25 


2.25 


2.99 


0.82 


138 




B 


2.30 


2.24 


2.37 


3.08 


0.88 


134 




A-A 


2.10 




2.24 


2.40 


0.24 


164 


1/4 


Ti-Ti 


2.15 


2.31 


2.30 


2.91 


0.55 


150 




T2-T2 


2.17 


2.28 


2.29 


2.92 


0.60 


148 




B-B 


2.20 


2.27 


2.37 


2.99 


0.71 


142 




Ag(110)p(2xl)O 


2.08 








0.10 


175 



10 



TABLE IV; Some key quantities of O/Ag(410) structures describing the bond lengths between Ag''^ and the nearest silver atom 
for various configurations (Ag^, Ag^ and Ag^ as in Tab. Illlll . The distances d^gX.^gY always refer to the distance between 
the two given silver atoms that are bonded to the same O adatom. The values in parentheses refer to the clean Ag(410) surface. 



e (ML) 


Configuration 


C^Aga-AgR (A) 


f^AgR-AgS (A) 


'^AgR-AgS' (A) 


C^AgR-AgB (A) 




A 


4.32 (4.16) 




3.08 (2.90) 


2.93 (2.84) 


1/8 


Ti 


4.26 


2.97 (2.91) 


2.98 (2.90) 


3.07 (2.91) 




Ta 


4.26 


2.98 (2.92) 


3.00 (2.92) 


3.04 (2.92) 




B 


4.24 


2.99 (2.92) 


2.98 (2.91) 


3.06 (2.96) 




A-A 


4.16 




3.15 


3.00 


1/4 


Ti-Ti 


4.16 


3.00 


2.99 


3.15 




T2-T2 


4.16 


2.99 


2.99 


3.12 




B-B 


4.16 


3.02 


2.98 


3.08 



TABLE V: Some key quantities describing the O chemisorption sites on Ag(210) in various geometries. The labels have the 
same meaning as in Tab. 11111 



e (ML) 


Configuration 


d-O-Ag^ (A) 


C^O-AgS (A) 


^-AgS' (A) 


do-AgB (A) 


Az (A) 


an 


1/4 


A 


2.21 




2.21 


2.42 


0.44 


158 




B 


2.31 


2.23 


2.40 


3.11 


0.92 


132 


1/2 


A-A 


2.09 




2.24 


2.42 


0.21 


168 




B-B 


2.26 


2.25 


2.41 


3.09 


0.88 


134 



TABLE VI: Some key quantities of O/Ag(210) structures describing the bond lengths between Ag''^ and the nearest silver 
atom for various configurations (Ag^, Ag^, Ag^ and Ag^ as in Tab. IllH . The distances d^gX-AgY always refer to the distance 
between the two given silver atoms that are bonded to the same O adatom. The values in parentheses refer to the clean Ag(210) 
surface. 



e (ML) 


Configuration 


C^AgR-AgR (A) 


f^AgR-AgS (A) 


'^AgR-AgS' (A) 


C^AgR-AgB (A) 


1/4 


A 


4.33 (4.16) 




3.04 (2.89) 


2.95 (2.85) 




B 


4.22 


2.96 (2.89) 


2.97 (2.91) 


3.04 (2.95) 


1/2 


A-A 


4.16 




3.09 


3.04 




B-B 


4.16 


2.95 


2.97 


3.03 




FIG. 1; Schematic representation of the Ag(410) (left) and Ag(210) surface (right). The labels A, B, Ti and T2 indicate 
different adsorption sites: A, above the step; B, below the step; Ti and T2, hollow sites at terrace. The surface unit-cell vectors, 
and the [001] and [InO], n = 4, 2, crystal axes are also shown. 
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TABLE VII: Vibrational frequencies and displacement eigenvectors for different O/Ag structures. The angles which define 
the directions of the eigenvectors are shown in Fig. [TTl For both the Ag(110)p(2xl)O and O/Ag(100)-c(2x2) systems 9 is the 
angle with respect to the axis perpendicular to the surface. For the Ag(110)p(2x 1)0 ip is the angle in the surface plane with 
respect to the added row. The values in parentheses are the LDA resultsi— 



Surface 


Configuration 


Mode 


V (meV) 


Q (°) 


<^ (°) 






TTll 


30 (35) 







AeCIOO) 


c(2x2) 


77I2 


50 (55) 


90 








7713 


50 (55) 


90 










28 


90 


90 




added row 


77I2 


38 











ma 


73 


90 









mi 


24 


23 


270 




A-A 


7TL2 


37 


70 


90 






77I3 


68 


90 





Ag(410) 




mi 


29 


9 


90 




T1-T1 


7712 


32 


81 


270 








58 


90 











24 


23 


270 




A"A 


m2 


37 


69 


90 






ma 


71 


90 









mi 


25 


30 


270 


Ag(210) 


A 


m2 


40 


64 


90 






ma 


46 


90 









mi 


27 


39 


90 




B-B 


m2 


34 


49 


270 








37 


90 
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Ti 
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Ti-Ti 
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To-B 
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FIG. 2: (Color online) Top views of optimized Si, Si-Si, and S1-S2 O/Ag(410) structures. Grey (red) balls represent silver 
(oxygen) atoms. 




FIG. 3: (Color online) Top view of the Ag(110)p(2x 1)0 surface. Grey (red) balls represent silver (oxygen) atoms. 
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FIG. 4: (Color online) A schematic definition of Ag^, Ag^, and Ag^ atoms, which are labelled as R, S (or S'), and B, respectively. 
Larger balls are Ag atoms, while smaller red ball is O adatom. The step-edge is marked with dot-dashed line. The Ag^ atoms 
are lying in the row of oxygen adatoms, while Ag^ atoms are the other nearest surface silver atoms. The Ag^ silver atom is 
just below the O adatom. 




FIG. 5: (Color online) Side views of the first three layers of Ag(410) substrate with O adatoms on site A in the (2x1) and 
(1x1) geometries. Grey (red) balls represent silver (oxygen) atoms. The inset in the top panel shows the side view of the 
optimized clean Ag(410). 



FIG. 6: (Color online) Same as in Fig. |S| but for T2 and T2-T2 configurations. 
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FIG. 7: (Color online) Side views of the Ag(210) substrate with O adatoms on site A (left panels) and on site B (right panels) 
in the (2x1) and (1x1) geometries. Grey (red) balls represent silver (oxygen) atoms. The inset in the top-left panel shows the 
side view of the optimized clean Ag(210). 




FIG. 8: (Color online) Charge difference density, i.e. the difference between the electron density of the O/Ag system and the 
density of the clean surface [A7i(r) = no/Ag('") — '^Ag('")], for A and A-A configurations (top views). Top panels: isosurfaces at 
—0.015 e/ao (blue). Bottom panels: charge density contours in a plane passing through O, Ag''' and Ag^ atoms. Contours are 
drawn in linear scale from —0.07 to 0.0 e/ap, with the increment of 0.01 e/ag. Grey (red) balls represent silver (oxygen) atoms. 
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FIG. 10: (Color online) ILDOS in the energy window (-7.0, -6.0)eV for the A-A configi 
at 0.015 e/ag (yellow) and the 10 contours are drawn in linear scale from 0.015 e/ajj to 0. 



juration (top view). The isosurface is 
.15 e/aj). 
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FIG. 12: (Color online) Displacement eigenvectors mi and m2 for the A-A and the Ti-Ti configurations on Ag(410). For the 
A-A (Ti-Ti) geometry the mode mi vibrates at 24 meV (29 meV) while the mode m2 at 37 meV (32 meV). 



